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Abstract. We translate the concept of succession rule and the ECO method into 
matrix notation, introducing the concept of a production matrix. This allows us 
to combine our method with other enumeration techniques using matrices, such as 
the method of Riordan matrices. Finally we treat the case of rational production 
matrices, i.e. those leading to rational generating functions. 
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1. Introduction 

In |DFRj a new type of infinite matrices, the so called production matrices, was 
introduced and studied. Such matrices are simply a new way to represent succession 
rules, which allows us to work with them using algebraic methods. 

Succession rules were first introduced by Julian West in [ WeltrWe2] . and later were 
used as a formal tool for the enumeration with the ECO method |BDLPP] (ECO 
stands for Enumeration of Combinatorial Objects). In a word, it is a constructive 
method to produce all the objects of a given class, according to the growth of a cer- 
tain parameter (the size) of the objects. Basically, the idea is to perform a prescribed 
set of operations (sometimes called local expansion) on each object of size n, thus 
constructing a set of objects of the successive size. This construction should induce 
a partition of all the objects of any given size (that is, all the objects of a given 
size are produced exactly once from the objects of immediately lower size, through 
the ECO construction). If an ECO construction is sufficiently regular, then it is 
often possible to describe it using a succession rule. Intimately related to the con- 
cept of succession rule is the notion of generating tree, which is the most common 
way of representing a succession rule. The principal applications of ECO method 
are enumeration |BDLPPT1 IBBDFG^ , random generation |BDLPP2] , or exhaustive 
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generation [BBGPl IDLFRj of various combinatorial structures. For all these topics 
we refer the reader to the rich survey |BDLPP] . 

In |DFR] we have proposed a translation of the concept of succession rule, and 
hence of the ECO method, into matrix notation. This is achieved by introducing 
two different possibly infinite matrices strictly related to a succession rule, namely its 
production matrix, and its ECO matrix. The main goal of our approach is to provide 
a representation of succession rules more suitable for computations; indeed we define 
some operations on production matrices in order to reproduce well known operations 
on the numerical sequences they represent. This leads to the determination of the 
generating functions of such sequences, often more easily than it was previously done 
by other methods (see |BDLPPl IBBDFGGl lFP] l 

In this paper we deepen the study of production matrices and ECO matrices, 
in particular trying to establish relations with some kinds of infinite matrices that 
have been recently introduced and studied, known as Riordan arrays |MRSVl [Rl IShl 
ISCWWl gg. 

In Section [3] a comparison with the Riordan matrices method is established, thus 
completing an investigation started in |MVj . In Section H] we outline the main results 
concerning exponential Riordan matrices |DSj . and then study this concept from the 
point of view of production matrices. In the last section we consider finite production 
matrices, representing finite succession rules. Throughout the whole paper, a huge 
amount of examples are described or simply sketched. 

Before ending this introduction, just a few words of explanation concerning some 
of the notations used in the work. We have selected a more or less standard symbol 
for the generating functions of frequently occurring sequences. Namely, 

m C = = 1 + z + 2z'^ + bz^ + Uz^ + 42z^ + 132^6 + 429^^ + ... for the 

Catalan numbers (A000108); 

• M = ^~^~^2~2^^~^^ - = 1 + 2 + 2^2 + + 9^4 ^ 21z^ + 51z^ + Vllz^ + ... for 
the Motzkin numbers (A001006); 

^ Yi= ^-^-^^-^z+z^ = i^2z + Qz'^ + 22^3 + 90^4 ^ 394^5 ^ 1306^6 + 8558^^ + . . . 

for the large Schroder numbers (A006318); 

• S = i+^-Vi-62+^^ = 1 + ^ + 3^2 ^ ^^^^3 ^ 45^4 ^ ^97^5 ^ 9Q3_26 ^ 4279^7 _^ _ _ 

42: 

for the small Schroder numbers (A001003); 

• T =1 + zT^ = 1 + z + ?>z^ + I2z^ + hhz^ + 273^^1428^^ + 7752z^ + ... for the 
ternary numbers (A001764); 

• F = /-^"lE^, = 1 + ^2 + 2^3 + Qz^ + 18^5 + 57^6 ^ ^gg^7 + ... for the Fine 

2(3— vl— 42) 

numbers (A000957); 
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• B = = l + 2z + 6z'^ + 20z^ + 70^^ + 252^^ + 924^^ + 3432z^ + ... for 

the central binomial coefficients (A000984). 

Throughout the paper the A****** number between parentheses following a se- 
quence is the identification number of that sequence in [Slj. Most of the matrices we 
are going to consider are infinite; their lines (rows and columns) will be indexed by 
nonnegative integers, and we will write "line 0" to mean the first line, "line 1" to 
mean the second line, and so on. 

For the sake of readability we have decided that the paper be self-contained. This 
implies, in particular, that some basic definitions and properties of production ma- 
trices, reported in Section 2, are taken directly from [DFRj . Therefore the reader 
familiar with these topics can skip Section 2. 

2. Basic definitions 

A succession rule is a formal system consisting of an axiom (a), a G N"*", and a set 
of productions: 

{(h) ^ ie^ih)){e2ih)) ■ ■ ■ (ekAh)) : t G N}, 

where Cj : N+ — > N+, which explains how to derive the successors 
(ei(/c)), (e2(A;)), . . . , (efc(A;)) of any given label (k), k G N+. In general, for a suc- 
cession rule Q, we use the more compact notation: 

^ ' i U {ei{k)){e,{k)) . . . {e,{k)) ' 

(a), (fc), (ej(A;)), are called the labels of Q (where a,k,ei{k) are positive integers). 
The rule fl can be represented by means of a generating tree, that is a rooted tree 
whose vertices are the labels of f2; (a) is the label of the root and each node labeled 
(k) has k sons labeled by ei{k), . . . ,ek{k) respectively, according to the production 
of (k) in ([T]). A succession rule fl defines a sequence of positive integers (/n)n>05 fn 
being the number of nodes at level n in the generating tree determined by Q. By 
convention the root is at level 0, so /o = 1. The function fn{x) = J2n>o fn^"" 
generating function determined by fl. 

In [DFRj we have proposed a representation of succession rules by means of infinite 
matrices P = {pk,i)k,i>o- Assume that the set of the labels of a succession rule is 
{{h)}k, and in particular that Iq is the label of the axiom. Then we define pk^i to be 
the number of labels k produced by label Ik- We call P the production matrix of the 
succession rule. Observe that the first row of a production matrix gives precisely the 
production of the axiom. 
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Labels do not occur explicitly in this matrix representation of succession rules. 
However, they are the row sums of P. In particular, the label /q of the axiom is the 
first row sum of P. 



Example. To the succession rule 

(2) 

{2k) (2)'=(4)---(2A;)(2A; + 2) 
there corresponds the production matrix 



(2) 



(3) 



P 



/I 1 

2 110 

3 1110 

4 11110 

5 11111 

V: ; ; ; ; ; 



(2) 




(2) (4) 




(2) (4) (2) (2) (4) (6) 




(2) (4) (2) (2) (4) (6)(2) (4) (2) (4)(2)(2) (4) (6) (2)(2)(2)(4) (6) (8) 

Figure 1. The first levels of the generating tree associated with the 
succession rule in 

In the generating tree at level zero we have only one node with label /o(= 2). This 
is represented by the row vector 

ro = (1 . . .) . 

At the next levels of the generating tree the distribution of the labels li,l2,--- is 
given by the row vectors Vi, i > 1, defined by the recurrence relation rj = ri^iP. 
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Stacking these row matrices, we 


obtain the matrix 




( 1 



















1 


1 
















3 


2 


1 











Ap = 


10 


6 


3 


1 










35 


20 


10 


4 


1 







126 


70 


35 


15 


5 


1 



V ; ; ; ; ; ; 

The row sums of the above matrix are 

1, 2, 6, 20, 70, 252, 924, 3432, 12870, 48620, . . . 

i.e. the central binomial coefficients. This is the sequence corresponding to the 
succession rule of our example. The enumerative properties of this succession rule 
have been examined in detail in |BFRj . Using the terminology of [DFR] . we will refer 
to as the ECO matrix induced by P. 

Remarks. Let P be the production matrix of a given succession rule VL. Through- 
out the whole paper we will denote 



u 



1,0,0,. ..,0,...) 



/1\ 

1 
1 



having appropriate sizes. The following properties have already been stated 
in [DFRj . 

(i) The labels of the nodes of the corresponding generating tree are the row 
sums of P. If two row sums happen to be equal, then, as labels, they will 
be considered to be distinct. This can be achieved by using, for example, 
distinguishing subscripts; in the vocabulary of succession rules, these are called 
colored succession rules (see [FPPR]). 

(ii) The distribution of the nodes having various labels at the various levels is 
given by the ECO matrix 

f \ 

u^P 



A, 



u^P 
\ ■■■ J 
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(indeed, we have tq = , ri = r^P = P. r2 = riP = P"^, . . .). The same 
fact can be expressed in a concise way by the matrix equahty 

(4) DAp = ApP, 

where D = (5i,j+i)i,j>o is the usual Kronecker delta). In some sources 
|PWt ISh] the matrix P is also called the Stieltjes transform matrix of Ap. 

(iii) The sequence a„ induced by the succession rule is given by a„ = u^P^e. 

(iv) The bivariate generating function of the matrix Ap is 



G{t,z) 



u 



\l-zP) 



-1 



/1\ 

t 

t^ 



VI 



The generating function of the sequence corresponding to the succession rule 
is 

fp{z)=u'{I-zP)-^e. 

The exponential generating function of the sequence corresponding to the 
succession rule is 



Fp{z) = exp{zP)e. 
Example. We intend to find the sequence determined by the production matrix 



(5) 



P 



/o 








Denoting by P„ the upper left n by n submatrix of P, it is not difficult to compute 
the exponential of the matrix zPn, since it is an upper triangular matrix. The eigen- 
values of Pn are easily seen to be the nonnegative numbers 0, 1,2, 3, — 1, each 
with multiplicity 1. Thus we have immediately 

exp{zPn) = Cexp{zDn)C~^, 

where D„ is the diagonalization of P„, so that exp{zDn) = {zSije'^^)o<ij<n, and C 
is a suitable invertible matrix. More precisely, a simple computation shows that C 
is an upper triangular matrix in which the entry has the form jp:jy^ (where 

obviously i < j). This implies that also is upper triangular and, for i < j, 
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its entry has the form (—1)-' *(jZi)T- Now the computation of the first row of 
exp{zPn) is immediate, and we find for it 



1 



i.(e^_l)2 l(e^-l)3 ... ^(e^-i)-i) 



Taking the sum of these entries and letting n ^ oo, for the exponential generating 
function induced by P we obtain 

Gp{z) = e'^-\ 

The corresponding sequence is 1,1,2,5,15,52,203,876,... (AOOOllO; Bell num- 
bers). 

3. RiORDAN PRODUCTION MATRICES 

An infinite lower triangular matrix A is called a Riordan matrix if its column k 
(fc = 0, 1, 2, ...) has generating function d{z)[zh[z)Y ^ where d{z) and h[z) are formal 
power series with (i(0) 7^ 0. If, in addition, /i(0) 7^ 0, then A is said to be a proper 
Riordan matrix. We may write A = {d{z), h{z)). 

Riordan matrices were first introduced in [SG WWj . where it has been proved that 
they constitute a natural way of describing several combinatorial situations. 

Proper Riordan matrices are characterized by a fundamental property found by 
Rogers in 1978 and then examined closely by Sprugnoli [Sp]. In fact, for a given 



proper Riordan matrix A = {d{z),h{z)) = {dn,k)n,k>o^ there exist unique sequences 
a = {ao, tti, a2, ■ ■ ■) (ao 7^ 0) and ( = (Co, Ci, C2) ■ ■ ■) such that 

(i) every element in column can be expressed as a linear combination of all 
the elements in the preceding row, the coefficients being the elements of the 
sequence (, i.e. 

(6) dn+lfl = Codnfl + Cldn,l + C2<^n,2 + • • • ; 

(ii) every element dn+i,k+i, not lying in column or row 0, can be expressed as 
a linear combination of the elements of the preceding row, starting from the 
preceding column on, the coefficients being the elements of the sequence a, 
i.e. 

(7) dn+l,k+l = OiQdn,k + OLidn,k+l + «2C^n,fe+2 + • • • 

Conversely, the existence of such sequences a and C, ensures that the matrix A is a 
proper Riordan matrix. 

The sequences a and C will be called the a-sequence and the ^-sequence of the 
Riordan matrix. By abuse of notation, by a, C we shall denote also the generating 
functions of these sequences. 
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The functions d{z), h{z), C{z), a{z) are connected by the relations 

h{z) = aizhiz)), diz) = ^-^j^^. 

It is also known [Sp] that, under the assumption d{0) = h{0) 7^ we have d{z) = 
h{z) if and only if a{z) = (io,o + ^Ci^)- 

At this stage, it is natural to investigate the relationship between the theory of 
production matrices and the theory of Riordan matrices. It turns out that, if the 
ECO matrix Ap, induced by a production matrix P, is Riordan, then the matrix P 
has a very simple structure. First we complete a result just outlined in Section [2l 

Lemma 3.1. Let P = {pi^k)i,k>o be an infinite production matrix and let Ap = 
{dn,k)n,k>o be its associated ECO matrix. Then, for any k,n > 0, we have that 
do,k = <^o,fc and 

(8) dn+l,k = dnflPo,k + dn,lPl,k + dn,2P2,k + dn,3P3,k + • • • , 

where 6 is the usual Kronecker delta. 

Conversely, if Ap and P are infinite matrices such that relations ^ hold, then Ap 
is the ECO matrix induced by P. 

Proof. The first part of the theorem is essentially statement (ii) in the Remarks of 
Section 2. The converse is easy. □ 

Proposition 3.1. Let P be an infinite production matrix and let Ap be the matrix 
induced by P. Then Ap is a Riordan matrix if and only if P is of the form 

...\ 
ao ... 

«! ao 
a2 cti ao 

OLZ OL2 OL\ 

a2 

■■.) 

Moreover, columns and 1 of the matrix P are the C- and a-sequences, respectively, 
of the Riordan matrix Ap. 

Proof. If P is as in ([2]), then from the relation r„ = r^-iP between the rows of the 
matrix Ap (see the above lemma) we can see at once that columns and 1 of P play 
the roles of the Q- and a-sequences, respectively. Consequently, by what has been 
said at the beginning of this section, the matrix is Riordan. Conversely, assume 
that y4p is a Riordan matrix. Then equalities ([6]) and ([7j) hold, and in view of the 



(9) 



/Co «o 

Ci «i 

C2 

Cs 

C4 

Cs 



a2 
as 

"5 
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above lemma they translate into the fact that the production matrix oi Ap must be 
as in ©. □ 

The above proposition, formulated in terms of succession rules, is the main result 
of [MY]. 

Because of the above property, a production matrix P having the form (Q will be 
called a Riordan production matrix. 

In the case of a given Riordan production matrix P, having (Cn)n>o and {an)n>Q as 
its first two columns, one can easily determine the bivariate generating function G{t, z) 
of the matrix Ap induced by P and then, obviously, also the generating function fp{z) 
of the sequence induced by P. 

Proposition 3.2. Let P he a Riordan production matrix and let C,{z) and a{z) he 
the generating functions of its first two columns, respectively. Then the hivariate 
generating function G{t, z) of the matrix Ap induced hy P and the generating function 
fp{z) of the sequence induced hy P are given hy 

d{z) 



(10) 



Gp{t,z) 



1 - tzh{z) 

where h{z) is determined from the equation 

(11) h{z) = a{zh{z)) 

and d{z) is given hy 

1 



d{z) 
- zh{z) 



12 



d{z) 



1 - zC{zh{z)) 



Proof. The proof is straightforward from general facts of Riordan matrices theory. 
The only thing which is worth remarking is that, in the case of ECO matrices, d{0) = 
1; this is the reason for which the numerator of the fraction in the r. h. s. of f|T2|) is 
1 (in general it should be d{0)). □ 

Example. Consider the Riordan production matrix 



/ 3 

7 

15 
31 
63 



1 
3 
7 
15 




1 
3 
7 



31 15 





1 
3 
7 









1 
3 1 



127 63 31 15 7 3 



v 
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Note that the row sums of P, i.e. the labels of the generating tree, are the Eulerian 
numbers 4, 11, 26, 57, 120, . . . (A000295). We have a{z) = (^^_,)l^_2,j and a - zC = 1. 
Then, recalling some of the results stated at the beginning of this section, is the 
Riordan matrix {d{z), h{z)) such that 

(13) h = — , d = h. 

^ ' {I - zh){l - 2zhy 

Let us denote by fp the generating function of the sequence determined by P. From 
the last theorem we get 

(14) 



h 



1- zh 

Eliminating d and h from equations f|T3l) . f|T^ . we obtain (1 + zfp)^ = fp{l — zfp). 
The substitution K = z + z^ fp leads to the equation = {K — z){2z — K), which is 
the equation giving the generating function for the number of noncrossing connected 
graphs. The sequence corresponding to fp starts 1, 4, 23, 156, 1162, . . . (A007297). 

Example. Consider the Riordan production matrix 



/ 


1 














2 


1 


1 











11 


3 


1 


1 








54 


12 


3 


1 


1 





272 


55 


12 


3 


1 


1 


1427 


273 


55 


12 


3 


1 



p 



where in columns 1,2,3,... we have the ternary numbers 2^^(^^) and in column 
the same numbers diminished by 1. Thus, a{z) = T{z) and C,{z) = — 
where the generating function of ternary numbers T{z) satisfies T{z) = 1 + zT^{z). 
Consequently, h = a{zh) = T{zh) = 1 + zhT^{zh) = 1 + zh"^, i.e. h{z) is the 
generating function of the sequence ■^^^{^n)'''^ ~ 0,1,2,.... After some elementary 
manipulations we obtain d = h{l — zh) and fp{z) = h{z). The sequence corresponding 
to fp starts 1,1,4,22,140,969,. . . (A002293). 

Remark. The last example can be easily generalized to construct a production 
matrix inducing the sequence j^^^z^y;^{^^) ,n = 0,1,2, .... 

In the next table we show numerous other examples. The first two entries in a 
row define a Riordan production matrix, the next two entries show the bivariate 
generating function of the induced matrix Ap and the generating function fp of the 
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1 

{l-z)(l-tz) 
l-tzT{z^) 

l-{l+t)z 
1 

{l-z){l-tzC) 
{l-z)(l-2z) 

(l-32 + z2)(l-2-t2) 
(l-32+22)(l-tz) 

M 
1-tzM 

F 
l~tzF 

(1-2)C^ 

1-tz 
C 

l-tzC 

(l-zF)C'^ 

1-tzF 
1+B 
2(l-t2) 
1 



1 

l-zT{z^) 

1-22 

C 

1-2 

1-2 
1-32+22 

1-2 
1-32 + 22 
M 

1-zM 

c 

C-1 
z 

C-1 

z 

C-1 

z 

1+B 



SC 



(jueiice 


A-iiuiiiber 


2,3,4,5,6,7,8,9,... 


AUUUUz / 


1,2,3,7,12,30,55,... 


A047749 


2,4,8,16,32,64,... 


A000079 


2,4,9,23,65,197,... 


A014137 


2,5,13,34,89,233,... 


A001519 


2,5,13,34,89,233,... 


A001519 


2,5.13,35,96,267,... 


A005773 


1,2,5,14,42,132,... 


A000108 


2,5,14,42,132,429,... 


A000108 


2,5,14,42,132,429,... 


A000108 


2,5,14,42,132,429,... 


A000108 


2,5,15,50,176,... 


A024718 


2,5,16,61,258,... 


A104858 


2,5,17,72,345,... 


A104859 


2,6,18,57,186,622,... 


A000957 


2,6,20,70,252,924,... 


A000984 


2,6,20,70,252,924,... 


A000984 


2,6,22,90,394,... 


A006318 


2,6,22,90,394,... 


A006318 


2,6,22,90,394,... 


A006318 


2,7,28,121,550,... 


A010683 


2,7,30,143,728,... 


A006013 


1,3,7,19,51,141,... 


A002426 


3,8,21,56,154,440,... 


A014318 


3,8,22,64,196,625,... 


A014138 


3,8,23,70,222,726,... 


A000782 


3,8,24,75,243,808,... 


A000958 


3,9,28,90,297,1001,... 


A000245 


3,10,35,126,462,... 


A001700 


3,10,35,126,462,... 


A001700 



1-2^ 
1 



1 

1-2 

(T3i)2 

1 + z 



1 
1 

1+2-22 
1-2 

1 + Z 

1 

l-2z 
1 

l-z 



1 + 2 
1-2 
1 

(1-2)2 

z 



2-2z+z2 

1-2 

2 

1-22 

2 

1- 2 

2- 2 
1-2 

2 

1-2 

2 + z 



1 

l + z 
1 

l-z 
l + Z 

1 

1 + z + 



2 


(1-2)' 


l-2z 


1-2^1 


C 
l-z 


1 


1 


1 


l-z 


1-2 


1 — 2 — 22 


(1-2)' 


g-2)(l-22) 


l-2z 

1 



1 + 2 
1-2 

(T^ 
1 

1-2 
1 

l-z 

(1 + ^) 

l-z 

1-22 

l+z 

1- z 

2- 2 

1-2 
1 + 2 



(l-2)(l-t2_R) 

1 

(l-2)(l-t2T2) 

FC 

1-tzC 
B 

C{l-tzC) 

l-te(C-l)) 
S 

1-tzS 
1+zR 

1-tzR 
1 

l-tz(l+R) 



l-z 

1 


(l-zR)(l- 
T 


-tzR) 


(1-2)2 


1-tzT'^ 




l + z + z"^ 

1 


M 




(1+2M)(1 

1 


-tzM) 


l-z 

1 


(1-22)(1- 

C 


tzC) 



l-z 

1 

l-z 

1 

l-z 

1 

l-z 

1 



{l-z){l-tzC) 
2C-1 



l-z 

{l + zf 



1 


-tzC 




F 


1 


-tzC 




C2 


1 


-tzC 




B 


1 


-tzC 




C2 


1 


-teC2 



2^-2) 
1-2 

T 
l-z 
F-1 

z'^ 

B 



B 
R 
R 
R 

g2 
rj-i2 



Vl-22-322 

C 

1-22 

_ci 

1-2 

2(g-l) 



2 

C-F 



BC 
BC 



-C 
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a{z) 



Gp{t,z) 





sequence 


A-number 


BC 


1,3,10,35,126,462,... 


A001700 


BC 


1,3,10,35,126,462,... 


A001700 


d2 

o 


1,3,11,42,163,638,... 


A032443 


5-1 

z 


1 3 11 45 1 Q7 Q03 


AOOl 003 


5-1 


1,3,11,45,197,903,... 


A001003 


T-1 

2 


1,3,12,55,273,1428,... 


AOOl 764 


T 


1,1,3,12,55,273,... 


AOOl 764 


C(22) 


1,3,13,67,381,2307,... 
1,3,13,67,381,2307,... 
1,4,15,56,210,792,... 


A064062 
A064062 
AOOl 791 


1-2(7(22) 

C(22) 


1-2C{22) 


C(32) 


1,4,25,190,1606,... 


A064063 


l-2C(3z) 



2-32-2^ 
(l-2)(l-22) 

1+C 

1-2 

2{l + z) 
1 

l-z 
2 

1-2 

2-2 
(1-2)2 

C-1 ]_ 

2 1 — 2 

2 

1-22 
1 

1-2 

3-22 

1-2 
1 

1-32 



1-22 



(1 + ^)2 

2-2 
1-2 
1 + 2 

1-2 
1 

(1-2)2 

c 

1 

1-22 

2 

1-2 
1 

1-2 

3 

1-32 



gC(l-2F) 
l-tzF 

{l-z)BC 
1-tz 
B 

l-tzC^ 

s 

l-tz{l+R) 

R 
l~tzR 

rp2 



l-tzT^ 
T(l~zT) 

l~tzT 
C(2z) 

l~tzC{2z) 



(l-zC(2z))(l-2tzC(2z)) 
BC 



l-tzC 
C(32) 

l-tzC(32) 



4. Exponential Riordan production matrices 

In this section we outline the main resuhs concerning exponential Riordan matrices 
as they are exposited in |DS] , and then study this concept from the point of view of 
production matrices. For the case of tridiagonal matrices, see also |All IA21 IPWj . 

Let d{z), h{z) be two formal power series such that d{0) 7^ 7^ h{0). An exponential 
Riordan (briefly, eR) matrix is an infinite lower triangular array A = {an^k)n,k>o whose 
column k {k = 0,1,2, .. .) has exponential generating function Ck{z) = yd{z){zh{z))'' . 
We will use the notation A = [d{z), h{z)] to denote the eR matrix determined by d{z) 
and h{z). From the definition it follows at once that the bivariate generating function 
of A, namely 

(15) GA{t,z) = y2an,kt''^, 

n,k 

is given by 

(16) GA{t,z) = d{z)exp(tzh{z)). 

Observe that GA_{t, z) has been defined to be ordinary with respect to the variable 
t and exponential with respect to the variable z. 

In |DS] it is shown that eR matrices form a group (with respect to the usual 
multiplication operation), whose identity element is [1, 1]. 



14 



EMERIC DEUTSCH, LUCA FERRARI, AND SIMONE RINALDI 



Looking at the above definition, one could expect many similarities with the theory 
of (classical) Riordan matrices. From the point of view of the present work, one of 
the main analogies with the ordinary case is the possibility of expressing every entry 
of an eR matrix as a linear combination of the elements of the preceding row. More 
precisely, we have the following result ^DSj. 

Proposition 4.1. Let A = {an,k)n,k>o = [d{z), h{z)] be an exponential Riordan matrix 
and let 

(17) c{y) = co + ciy + C2y'^ + . . . , r{y) = tq + riy + r2y'^ + . . . 

be two formal power series such that 

r{zh{z)) = {zh{z)y 
<^h{z)) = ?g ■ 

Then 

(19) (i) an+1,0 = y^j^-Cian,i 




(20) 



0'n+l,k — r oa„ 



i>k 



kri. 



-k+ljO-r, 



or, defining C-i = 0, 



(21) 



■ i>k-l 



+ kri 



i-k+l 



Conversely, starting from the sequences defined by [T^ , the infinite array {an,k)n,k>o 
defined by l[21\) is an exponential Riordan matrix. 

Remark. The sequences {cn)n>o, (?'n)n>o are called respectively the c-sequence and 
the r-sequence of A. There is a clear analogy with the a-sequence and the C-sequence 
of the classical case. However, whereas for an ordinary Riordan matrix the coeffi- 
cients of the linear combinations in ([7]) do not depend on the column index, in the 
exponential case they do, as it is clear from formula (pTj) . 
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Consider the infinite matrix P of tlie coefficients in ( 1211) 

..A 




ro 

j(co + 4ri) ... 

•••/ 

or, in a compact form, 

i\ 

(23) P = {Pi,j)i,j>o, Pi,j = —{ci^j + jvi^j+i) (c_i = 0), 

so tliat formula ( l2Ti) can be rewritten as a matrix equality: 

AP = DA, 

where, as in formula (Hj), D = (5ij+i)jj>o {S is the usual Kronecker delta). The 
matrix P is "almost" lower triangular, meaning that Pi+i^i = tq; whereas Pi+k,i = 
for > 2. 

Denoting by diag(n), n > —1, the sequence {pn+k,k)k>o, the matrix P is character- 
ized by the following facts: 

• diag(-l) is constant ro, ro, ro, . . ., 

• diag(O) is an arithmetic progression with first term cq and ratio ri, 

• diag(l), after division by 1,2,3,..., is an arithmetic progression with first 
term Ci and ratio r2, 

• diag(2), after division by 1 ■ 2, 2 ■ 3, 3 ■ 4, . . ., is an arithmetic progression with 
first term C2 and ratio r^, etc. 

If the matrix P has nonnegative integer entries, then it can be viewed as a produc- 
tion matrix. In this case it will be called an exponential Riordan production matrix. 
Clearly, in order that [d{z), h{z)] be an ECO matrix, we have to set d{0) = 1 as the 
remaining initial condition of the system f|T8|) . 

It is convenient to have an expression for the bivariate exponential generating 
function of the matrix P, namely 

^-^ nl 

n,k 

Observe that again ipp has been defined to be ordinary with respect to the variable 
t and exponential with respect to the variable z. By simply replacing the values of 



(22) 



P 



/Co ro 

l!ci YT(co + ri) ro 

2!c2 fT(ci + r2) |(co + 2ri) 

3!c3 ft(c2 + r3) |(ci + 2r2) 

4!c4 f|(c3 + r4) |(c2 + 2r3) 






||(co + 3ri) 
|(ci + 3r2) 
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the Pn,kS in ( HM we immediately obtain: 



n.k 



n,k 



J2 -fcT'^"-'^"'' + ^ 



n,k 



Z' 

n\ 
^k^k 



n—k 



k^k 



^k^k 



k\ 



k\ 



r„z 



I.e. 



(25) 



ipp{t,z) = e'\ciz) + tr{z)). 



Clearly, setting t = 1 gives the exponential generating function of the row sums of 
P. If P is an eR production matrix, these are the labels of the succession rule induced 
by P, and their exponential generating function is: 

ipp{l,z) = e'{c{z)+r{z)). 

We consider in detail two examples. In the first one we are given an eR matrix 
A and we will look for a production matrix that induces A, while in the second 
example we start with an exponential Riordan production matrix P and we will find 
the induced ECO matrix. 



Example. Consider the infinite matrix of the unsigned Stirling numbers of the first 
kind 



A 



\ ; ; ; ; ; ; •./ 

Its bivariate exponential generating function is exp(— t log(l — z)) and, con- 
sequently, d{z) = and h(z) = — iHiiizfl. Now from the system (fT8l) we obtain 
r(— log(l — z)) = c(— log(l — 2;)) = leading to r{y) = c{y) = exp{y). In order 
to identify the corresponding matrix P, we find from fl25|) its bivariate exponential 
generating function ipp{t, z) = {1 + t) exp((l + t)z). Now it follows at once that P is 



/ 1 

















1 


1 














2 


3 


1 











6 


11 


6 


1 








24 


50 


35 


10 


1 





120 


274 


225 


85 


15 


1 
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the matrix of the binomial coefficients: 



P 



/I 

1 

1 
1 
1 
1 



1 
2 
3 
4 
5 





1 
3 1 
6 4 
10 10 









1 

5 




1 



6 15 20 15 6 



and the labels are 2, 4, 8, 16, 32, . . .. 



Example. Consider the matrix of the falling factorials (permutation coefficients): 



/ 1 


1 














2 


2 


1 











6 


6 


3 


1 








24 


24 


12 


4 


1 





120 


120 


60 


20 


5 


1 



The labels (i.e. the row sums of P) represent the total number of arrangements of 
a set with n elements (A000522). 

Looking at the diagonals of P, it follows at once that P is an eR production matrix; 
namely, the one determined by 



System (|T8l) becomes 

{zh{z)y = 



r{y) 



1 



1 



d'iz) 



1 



i-zh{zy d{z) {i-zh{z)y' 

with initial conditions h[0) = 1, d{0) = 1 and we obtain 

1 - VI - 2z 



d{z) 



v/r^2i ' 



h{z) 
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SO that Ap 
to be 



1 



Using the definition of an eR matrix, Ap turns out 



2 



Ap 



/I 

1 1 

3 3 1 

15 15 6 10 

105 105 45 10 1 

945 945 420 105 15 1 



\ 



The rows of are the coefficients of the Bessel polynomials with exponents in 
decreasing order (see A001497; the sequence of the row sums of Ap is A001515). 

In closing this section we remark that it is easy to pass from an ordinary Riordan 
matrix {d{z),h{z)) to the exponential Riordan matrix defined by the same pair of 
formal power series, that is [d{z), h{z)]. Indeed, if {d{z),h{z)) = {an,k)n,k>o and 
[d{z),h{z)] = (a„,fc)n,fc>o, then we have 



A succession rule is said to be finite if it has a finite number of productions. This 
property simplifies the structure of the rule, and many problems still open for the 
whole class of succession rules turn out to be solvable for the class of finite succession 
rules. In |BDPR] it is proved that the equivalence problem, i.e. the problem of estab- 
lishing if two given rules define the same numerical sequence, is decidable for finite 
succession rules. Moreover finite succession rules lead to rational generating func- 
tions, and have a combinatorial interpretation in terms of regular languages [DFPRj . 
More specifically, defining 

(1) TZ the set of rational generating functions of integer sequences (Z-rational 
functions, in the notation of [SSJ); 

(2) 7^+ the set of rational generating functions of positive integer sequences; 

(3) REG the set of generating functions of regular languages; 

(4) S the set of rational generating functions of succession rules; 

(5) JF the set of generating functions of finite succession rules, 

from [SSj and |FPPRj we obtain the following inclusions: 




5. Finite and rational production matrices 
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REG 

s 

The classes 71, REG, and are decidable, while 7l~^ is not decidable. It is also 
conjectured that J-' = S, i.e. every rational succession rule (which is a succession rule 
having a rational g.f.) is equivalent to a finite one. 

One of the most popular finite succession rules is the following, defining Fibonacci 
numbers, 1,1,2, 3, 5, 8, 13, 21, . . . (sequence A000045): 

;i) 

(26) { (l)-(2) 

(2)-(l)(2). 

A production matrix associated with a finite succession rule has finite dimensions: 
the number of rows is equal to the number of productions of the rule, whereas the 
number of columns is equal to the number of labels. For example, the production 
matrix associated with the rule in (!26l) is 

^0 r 

We recall that this example was also considered in [Welj . 

In this section we will determine some properties of finite and rational production 
matrices, re-obtaining and sometimes extending some results in |BBDFGG] . |BDPR] . 
and [FPPRj . 

Proposition 5.1. If P is a finite production matrix, then the induced generating 
function fp{z) is rational. 

Proof. This follows at once from fp{z) = u^{I — zP)~^e. 

Proposition 5.2. Let P be a finite production matrix and let 

+ Clt''-^ + ...+ Cg-lt + Cg 

be the minimal polynomial of P. Then the sequence (a„) induced by P satisfies the 
linear homogeneous recurrence equation 

an + Cia„_i + ... + Cg_ia„_g+i + Cqan-q = 0. 
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Proof. The matrix P satisfies 

pi + dP'^-i + ... + Cg^lP + Cgl = 0, 

from wliere, premultiplying by P'^''^ and postmultiplying by e, we obtain tlie con- 
clusion of tfie theorem. 

Remarks. 1. Clearly, this last theorem implies again that the generating function 
induced by a finite production matrix is rational. 

2. The last theorem can be slightly improved. Instead of the minimal polynomial 
one can take the P-annihilator of e, i.e. the unique monic polynomial which generates 
the ideal of all polynomials g{t) such that g{P)e = 0. This is a divisor of the minimal 
polynomial of the matrix P (see, for example, [HKj . p. 228). However, even this need 
not be the recurrence equation of lowest order satisfied by the sequence induced by 
P. 

Example. Consider the succession rule defined by the production matrix 



P = 



(2 


1 


1 


o\ 





3 











1 


2 


1 




1 


1 


V 



The minimal polynomial of P is t'^ — lOt'^ + 36t^ — 55t + 30, yielding the recurrence 
equation of order four: 

an — lOort-i + 36a„_2 — 55a„_3 + 30a„_4 = 0. 

However, the P-annihilator of e is t'^— 8t^+20t— 15, yielding the recurrence equation 
of order three: 

a„ - 8a„_i + 20a„_2 - ISa^-s = 0. 
This equation is satisfied by the sequence of the z-th row sums of the matrices 
J, P, P^, where i = 1,2, 3, 4. But, in this case, the sequence of the first row sums 
of the matrices I,P,P^, i. e. the sequence induced by P, satisfies the recurrence 
equation of order two: 

a„ - 5a„_i + 5a„_2 = 0. 
(the polynomial — 5t + 5 is a divisor of — 8t^ + 20t — 15). 

In the last part of this section we will focus on rational production matrices, i.e. 
production matrices whose induced generating function is rational. The conjectured 
equivalence between finite and rational succession rules is the motivation for the 
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present investigation. We will just state a simple result on rational production ma- 
trices along with some examples. The rational production matrices considered in 
such examples are not finite; in some cases the determination of an equivalent finite 
production matrix is a challenging problem. 

Proposition 5.3. Let P be a production matrix. If the vectors e, Pe, P^e, P^e, ... are 
linearly dependent, then the induced generating function fp{z) is rational. 

Proof. The hypothesis says that 

aoe + aiPe + ... + auP^e = 
for some A; G N and ao, ...,ak € R. Premultiplying by u^P" we get 



for any n e N. Therefore (an)n>o satisfies a linear recurrence, which implies that 
its generating function is rational. □ 

Remark. The above theorem generalizes Proposition 2 of [BBDFGG]. Indeed, the 
linear dependence P^e = aPe + (3e is equivalent to the "affine function" assumption 
of Proposition 2 of [BBDFGG]. 

Example. We consider the succession rule 

[(2){(A;) ^ (2)^-2(2 + niod 2)){k + 1)}] 
(see [BBDFGG]). Its production matrix is 



Introducing the vectors 



w — 



(0 

(1 





A 


1 











.. 






1 


1 


1 








.. 






3 








1 





.. 






3 


1 








1 


.. 






5 














1 . . 






5 


1 











.. 






V 












7 


1 





1 





1 





1 


1 


2 


3 


4 


5 


6 


7 


8 9 


10 



we can easily derive 



Pe — e + w, Pv — e, Pw — e + v + 2w, 
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from where 

Pe = e + w, P'^e = 2e + v + 3w, P^e = 6e + 3f + Sw. 

Eliminating v and w from the last three equalities, we obtain 

P^e - 3P^e + Pe - e = 0. 

Consequently, the generating function induced by the matrix P is rational. This 
generating function can be found by the standard procedure from the recurrence 
equation that follows from the previous relation, namely a„ — 3a„_i + o„_2 — a„_3 = 
and the initial values ao = l,ai = 2, 02 = 5 (the first components of the vectors 
e, Pe, P^e, respectively). One finds 

•^^(^) ^ 1 — o — r~2 — 3- 

1 — 3^ + 2:^ — 2;'* 

It is easy to observe that the rational production matrix P is equivalent to the 
following finite production matrix: 




Remark. In an entirely similar manner we can derive that the succession rule 
[(2){(k) ^ (2f-\3-(k mod 2))(k + l)}] 
(see [BBDFGG]) induces the generating function YZ^zf^lf^- 

Example. We consider the succession rule defined by the production matrix 

/I 1 

3 1 

6 1 

10 1 

15 1 

21 



V 



where in the first column the i-th entry is ^^^ji^. Denoting by [q] a column vector 
with i-th entry Cj, we have 

P[l] 











3. " 






5. " 








+ -. + 1 








2 






2 






2 


+ -. + l_ 
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By linearity we obtain 

Pe = l\^'] + l\^]+e, P^e = ^[.^] + ^[.] + 2e, P^e = ^[.^] + y [^] + 6e. 

Eliminating [i] and [P] from the last three relations, we obtain 

P^e - AP^e + 3Pe - e = 0, 
leading to the recurrence equation 

a„ - 4a„_i + 3a„_2 - ctn-s = 0. 

This equation with the initial values Oq = l,ai = 2,a2 = 6, obtained from the first 
entries of e, Pe, and P^e, respectively, yield the generating function 

f (z) = ^^-^)' 

■' ^ ' l-4:Z + 3z^ - z^' 

The first ten terms of the sequence are 1,2,6,19,60,189,595,1873,5896,18560... (se- 
quence A052544 in H]). 

Remark. If in the previous example, in the first column of P we take, more generally, 
the i-th term to be ai'^ + Pi + j, then, in the same manner, for the recurrence equation 
we obtain 



a„ - (a + /? + 7 + 3)a„_i - {a - (3 - 2-f - 3)a„_2 - (7 + l)an-j 



0. 



Example. We consider the succession rule defined by the production matrix 



P 



a 


1 














a 


1 


1 











a + 1 








1 








a + 1 


1 








1 





a + 2 














1 


a + 2 


1 















v 



where a is a nonnegative integer. It is easy to see that 

Pe= {a + l)e + v, 

where v = (0 112233445 ...)^. We also have 

Pv = e + V. 

Eliminating v from the last two relations, we obtain 

P^e - (a + 2)Pe + ae = 0, 
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leading to the recurrence equation 

a„ — (a + 2)a„_i + aan-2 = 0. 

Making use of the initial conditions Oq = l,ai = a + 1, we obtain the generating 
function 



^ ^ 1 - {a + 2)z + az^' 

For a = we obtain the sequence 1, 1, 2, 4, 8, 16, ... (A000079), for a = 1 we obtain the 
odd-indexed Fibonacci numbers 1, 2, 5, 13, 34, 89, 233, ... (A001519), and for a = 2, 3, 4 
we obtain sequences A007052, A018902, A018903, respectively. 
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